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We study the topological properties of full QCD with four flavours of dynamical staggered fermions. In particular the 
topological susceptibility is measured and the problem of the determination of its first derivative is discussed. 



1. INTRODUCTION 

Topology plays an essential role in determining the 
low energy features of strong interaction physics. In 
this paper we present a lattice study of the topologi- 
cal properties of QCD in presence of four flavours of 
degenerate dynamical staggered fermions. 

A relevant quantity is the topological susceptibility, 
defined in the continuum as 
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where K^(x) is the Chern current and Q(x) 
d^Kfj,(x) is the topological charge density, 



Q(x) = 
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Eq. Q defines the prescription for the singularity of 
the time ordered product when x — > . 

The value of x m tne quenched theory is related 
to the rj' mass by the Witten-Veneziano formula. It 
has been successfully measured on the lattice 1 2 1, con- 
firming the Witten-Veneziano prediction. 

In full QCD with spontaneous chiral symmetry 
breaking, \ is related to the quark condensate [ 3 1 



=0 + o(m q ) 



(3) 



where m q is the quark mass and Nf the number of 
flavours. We want to test Eq. (3) for Nf = 4 and 
different values of m a . 
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Another relevant quantity is the slope at q 2 = of 
the topological susceptibility, x', defined as 
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where x(<7 2 ) = / d 4 x e^ x (Q(x)Q(O)). In the 
quenched theory the consistency of the Witten- 
Veneziano mechanism requires \' to t> e small. In- 
stead in full QCD it is expected to be larger and re- 
lated to the singlet axial charge (and so to the proton 
spin crisis) |4|. Therefore its determination is of par- 
ticular importance. Unfortunately, techniques which 
work well for the lattice determination of x> are not 
straightforwardly applicable to the measurement of 
x' ■ We will present a preliminary estimate of x' and 
discuss perspectives for future refinements. 

2. TOPOLOGICAL SUSCEPTIBILITY 

We have simulated the full theory using the HMC 
algorithm with 4 flavours of staggered fermions at 
(3 — 5.35 and four different values of the bare quark 
mass, a • TYiq — 0.010, 0.015, 0.020, 0.050, perform- 
ing for each mass value respectively 6000,1500,1000 
and 3000 units of molecular dynamics time. We have 
used a 16 3 x 24 lattice and the Wilson action for the 
pure gauge sector. 

The field theoretical method has been used to deter- 
mine x- Given a discretization Q l {x) of the topologi- 
cal charge density, we define the lattice susceptibility: 
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Ql is related by a finite multiplicative renormaliza- 



Figure 1 . Topological susceptibility in lattice spacing 
units at different bare quark masses and j3 — 5.35. 

tion Z to the continuum topological charg^. More- 
over xl in general does not meet the continuum pre- 
scription used in Eq. (1) and, due to contact terms in 
the product Q(x)Q(0) at small x, a further additive 
renormalization appears 



XL = a 4 Z 2 ([3)x + M({3) . 



(6) 



The heating method |6| has been used to evaluate 
the renormalizations. Z is determined by measuring 
(Ql) on a sample of thermalized configurations ob- 
tained by heating a semiclassical one-instanton con- 
figuration. Similarly M is obtained by measuring xl 
on a sample with zero topological background, using 
the fact that Eq. (1) leads to x — in the trivial topo- 
logical sector. Smearing techniques have been used 
to improve the operator Ql 0> in order to reduce the 
renormalizations and obtain a better accuracy on x af- 
ter subtraction. We have performed 2 smearing steps, 
starting from a standard discretization of Q, with a 
smearing coefficient c = 0.9, as reported in 1 2 1. 

In Fig. 1 we report the results obtained for a 4 x ver- 
sus the quark mass. The clear dependence of a 4 x on 
m q disappears when going to physical units. We have 
fixed the scale by measuring the string tension, ob- 
taining a 2 (j = 0.073(7) at a ■ m q = 0.05 and a 2 a = 
0.033(4) at a ■ m q = 0.01. From these values, assum- 
ing \fo = 440 MeV, we obtain a(a ■ m q = 0.01) = 
0.081(5)fm and a(a ■ m q = 0.05) = 0.0121(6)fm, 
from which (x) 1/4 (a-™ g = 0.01) = (153±16) MeV 
and (x) 1/4 (a • m q = 0.05) = (153 ± 11) MeV. 



3 In the full theory mixings of Ql with fermionic operators ap- 
pears, which anyway can be shown to be negligible (5). 



We cannot rule out possible systematic errors com- 
ing both from the poor sampling of topological modes 
at the lowest quark masses |8| and from the deter- 
mination of the physical scale. Indeed at a ■ m q — 
0.01 we have determined the lattice spacing also 
by measuring m p and m % , obtaining a different re- 
sult, a — 0.101(5)fm. Using this value we obtain 
(x) 1/4 (a • m q = 0.01) = (123 ± 10) MeV, which is 
closer to the theoretical expectation coming from Eq. 
(4), (x) 1/4 ^ 110 MeV. 

3. ESTIMATE OF x ' 

On the lattice we can define 

.,2 



xi = ~£(<M<Mo)>z 2 
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A relation similar to Eq. (6) holds in this case |9 |: 

X ^a 4 Z 2 (/3) X ' + M'(/3). (8) 

While Z in the previous equation is the same ap- 
pearing in Eq. (6), since it is purely related to the 
renormalization of the topological charge, M' is a 
new additive renormalization, containing mixings of 
Xl to operators of equal or lower dimension. More 
generally, defining the two-point correlation func- 
tion of the topological charge density on the lattice, 
(Ql{x)Ql(0)}, one can write 



(Ql(x)Ql(0)) = a*Z 2 (Q{x)Q(0)) + m(x) 



(9) 



m(x) indicates mixings with terms of equal or lower 
dimension in the Wilson OPE of (Ql(x)Ql(0)), an d 
is related to M and M' in the following way: 



M = 
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M' = rn{x)x 2 



(10) 



M' cannot be determined by the heating method, 
since in this case the continuum x' is n °t constrained 
to be zero in the trivial topological sector. Therefore 
the techniques used for x cannot be straightforwardly 
applied to the determination of x' ■ I n principle M' 
can also be computed by lattice perturbation theory, 
but in practice this approach is not particularly suc- 
cessful, especially when dealing with smeared opera- 
tors for which the convergence properties of the per- 
turbation series worsen. We stress that also cooling 
techniques, which can usually be used to determine x> 
fail in the determination of x'- only the zero-moment 
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Figure 2. (Ql(x)Ql(O)) for the operators after one 
and two smearing steps, at a ■ m q = 0.01 

of the two-point function (Q(x)Q(O)), i.e. is topo- 
logically protected, since it can be expressed in terms 
of the global topological charge, \ — (Q 2 ) /V, an d Q 
is quasi-stable under cooling. This is not the case for 
higher moments, and in particular for x' ■ 

One possibility is to use an improved topologi- 
cal charge operator for which M ', even if not com- 
putable, is known to be small and thus negligible in 
Eq. (8). If we follow this ansatz for the 2-smeared 
operator used previously we obtain, from our simula- 
tion at a ■ m q = 0.01, ylxl ~ 20 MeV, in good 
agreement with the value expected from sum rules 
S. \f\x f \ = (25 ± 3) MeV. However, the system- 
atic error deriving from neglecting M', M'/x', can 
be estimated to be of the same order of magnitude as 
M/x — 40% at a ■ m q = 0.01, that is still quite large. 

A better estimate of x' requires deeper knowledge 
of the two-point correlation function. In the con- 
tinuum (Q(x)Q(0)) is known to be negative, by re- 
flection positivity, for x > 0, and positive and sin- 
gular for x — 0. Similarly, when using a lattice 
action which preserves reflection positivity, we ex- 
pect (Ql{x)Ql(0)) < whenever the two opera- 
tors Ql(0) and Ql(x) do not overlap. This is clear 
in Fig. 2, where a determination at a ■ m q — 0.01 
°f (Ql(x)Ql(O)) averaged over spherical shells of 
width Sx = 0.6a is reported for 1 and 2 smearings. 

The information on (Ql(x)Ql(O)) is not enough 
to extract (Q(x)Q(O)): according to Eq. (9), also 
m(x) is needed. However we know that m(x) comes 
from contact terms, so it must be zero at large |x|. 
Therefore, for a given operator Ql(x), there must be 
an xq such that, for |x| > Xq, m(x) can be ignored in 



> 
o 



1 smearing 
▼ 2 smearing 



0.2 0.3 0.4 0.5 0.6 0.7 

Ixl (fm) 

Figure 3. — (Q(x)Q(Q)) from 1 and 2-smeared oper- 
ators at a ■ m q = 0.01 

Eq. (9) and (Q(x)Q(O)) can be easily extracted, the 
value of Z being known from the determination of x- 
A practical way to extract Xq is the following: 
(Ql(x)Ql(0)) is determined for two operators Qli 
and Ql2 and one looks for a plateau at large |x| 
in the ratio of the two functions, corresponding to 
the squared ratio of the multiplicative renormaliza- 
tions, (Z1/Z2) 2 . In this way we have determined 
(Q(x)Q(0)) for 1 and 2 smearing steps at large |x| 
and a ■ m q = 0.01, as shown in Fig. 3. There is 
good agreement between the two determinations, as 
expected. Work is in progress to extrapolate the infor- 
mation at large |x| to smaller distances, thus allowing 
a more careful determination of x' ■ 
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